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Abstract. The results of the 2D modelling of Supernova 
remnant (SNR) RCW86 are given. Models of this rem- 
nant, which for the first time interpret the anisotropy of 
surface brightness as a result of the evolution of adiabatic 
SNR in the interstellar medium with a large-scale gradient 
of density, are considered. Estimations on the basic char- 
acteristics of RCW86 and the surrounding medium which 
follow from these models are found. It is shown that the 
observed surface brightness distribution of RCW86 may 
be obtaned in the both proposed up till now models with 
different initial assumptions: one about a Supernova ex- 
plosion in 185 A.D. and another about an explosion in the 
OB-asossiation. In order to obtain the observational con- 
trast of surface brightness it is necessary to have a medium 
with the characteristic scale of nonuniformity 11 pc if the 
age of RCW86 is 1800 years or 20 - 25 pc when the SNR 
is distant from us on 2.8 kpc. The preshock density con- 
trast between the southwestern and northeastern parts of 
RCW86 is in range 3.5 - 4.5. 

Key words: ISM: supernova remnants - ISM: individual 
objects: RCW86 - hydrodynamics - X-rays: interstellar 



1. Introduction 

Supernova remnant (SNR) G315.4-2.3 (RCW86, MSH14- 
65) is the result of a supernova explosion which is regarded 
as one of the most possible candidates for the " guest star" 
of 185 A.D. acco rding to Chinese manuscripts (Clark & 
Stephenson 1977). 



Chin & Huang ( |1994| ), Schaefer ( |l995D argue that the 
"guest star" in 185 A.D. was not a Supernova (SN). Then 
it is not right to assume that the age of RCW86 is t — 



1800 years although the SNR might be a relatively young 
remnant because it has a radial oriented magnetic field 



1994) 



(Milne 1987) hke Tycho, Cas A and Kepler SNRs (Strom 



The SN progenitor might be a member of OB- 
assosiation. In this case we may expect the SN to have 
been of type II, but it could not be born in 185 A.D. since 
the SNR has much bigger size than the expected one if it 
is situated at the distance d « 2.8 kpc obtained from the 
kinematic survey of ionized H (Rosado et al. 1996| ). 

From the ratio Ha/Hfj for optical filaments 



Ruiz 
1 kpc. 



( 1981 ) estimates the distance to the SNR as d ' 

R CW8 6, as radio source MSH14-65', was identified by 
Hill (1967). E — Z? relation yields an estimation on the 
distance to the SNR d = 2 3.2 kpc (Milne |1970D , but 
the possible range is, nevertheless, ci = 1 10 kpc (Strom 
199^. 



Both the radio and the X-ray observations reveal that 
the shape of RCW86 is close to a spherical one with the 
size 43' X 40' (at frequency 843 MHz) and is of a shell- 
like type with the approximately axis-symmetrical surface 
brightness distribution and higher emission in the south- 
we st (SW ) part of the SNR (CasweU et al. |l975i Pisar ski et 
al. |1984[ Claas et al. |l989l Whiteoak & Green |l996[ ). The 
brightest optical filaments coincide with the maximums of 



emission in radio and X-ray bands (van den Bergh 1978 
Pisarski et al 



1984) 



The important results in investigation of RCW86 were 
caused by analyzing the spectral properties of X-ray emis- 
sion from this SNR. The first soft X-ray observations 



were r eported and interpreted by Naranan ( 1977 ) . Winkler 
(1978), having used hard X-ray data, showed that the ob- 
served spectrum might be explained as thermal, with two 
bremsstrahlung components (from the forward and reverse 
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shock waves) with the temperatures Tiow ~ 0.22 keV and 
Thigh ^ 5 keV. 

The first map of the X-ray surface brightness distri- 
bution of RCW86 was presented by Pisarski et al. (1984) 
(EINSTEIN observatory). The authors noted that higher 
emission from the SW part might be caused by the in- 
teraction of the shock wave with the interstellar medium 
(ISM) with a density 2-3 times higher than the aver- 
age. The one temperature fitting of the spectrum from 
the whole remnant (it was accepted that the plasma is 
uniform and in collision ionization equilibrium (CIE)) 
yielded temperature T = 1.2 keV and column density 
Nu = 2.8-102" cm-2_ ^he obtained parameters of RCW86 



into account both the NEI plasma emission effects and 
the ISM density gradient. Therefore, previous studies of 



RCW86 are based on the one-dimentional Sedov (1959) 



are given in Table |^. Pisarski et al. (1984) also showed 



that if this SNR is at the adiabatic stage of its evolution, 
it cannot be the result of the type II SN. 

Nugent et al. ( |1984D analysed the spectrum of RCW86 
(HEAO 1 experiment) with the assumptions: 1) that 
plasma emits in the CIE and, 2) for the first time, that 
the plasma is under nonequilibrium ionisation condition 
(NEI) with Te = T,. The best fitting of the CIE spectrum 
is for Thigh 5.1 ± 0.14 keV, Tjow = 0.52 ± 0.04 keV, 
Nu = (1.1 ± 0.3) • 10^1 cm^2_ The NEI model gives 
A^H = (4.4 ± 0.3) • 10^^ cm~2 and the parameters of the 
SNR shown in Table |. The distance to RCW86 is smaller 
than to the OB-assosiation but the authors have noted 
that the possible deviation from the NEI model may in- 
crease the distance to d = 2.5 kpc. 



Claas et al. ( |1989D also used both the CIE and the NEI 
model for the interpretation of EXOSAT observation of 
RCW86. The possible contamination by the galactic ridge 
which increases the flux above 6 keV was subtracted, so 
in the two-component spectral model the smaller Thigh = 
3.4 ± 0.2 keV was obtained. The low temperature com- 
ponent is believed to be a NEI effect. In this paper the 
parameters of RCW86 were also estimated using the sur- 
face brightness distribution. Higher emission from the SW 
part m ust be caused by a more dense ISM. Therefore, the 
Sedov ( 195£ ) solution is used only for the interpretation 
of the northeastern (NE) part of the SNR. The NEI model 
yields Ts = Te/1.3 = Thigh/1.3 = (3.03 ± 0.2) • 10^ K. 

ASCA X-ray data of RCW86 are considered by Vink et 
al. ( 1997 ). The authors reveal a remarkable temperature 
variation over the SNR: T^ = 0.8 5 keV. They estimate 
the preshock density contrast for the SW and NE parts of 
RCW86 as 5 ^ 50 times. 

Modelling the SNR emission it is necessary to take 
into consideration the possibility that the plasma may be 
under the NEI conditions, since the effects of the NEI are 
important. At the same time the nonuniformity of the sur- 
rounding ISM is very important, too. Calculations show 
that the influence of a nonuniform medium on the evolu- 
tion and the properties of the emission of adiabatic SNRs 
may exceed the in fluenc e of the NEI in 10^ 10^ times 
(Hnatyk & Petruk 199^ , hereafter Paper I). At present it 
is practically inpossible to construct models which take 



solution which cannot restore the observed morphology of 
the SNR. In this paper we show that the global anisotropy 
of the surface brightness distribution of RCW86 may be 
explain as a consequense of the SNR evolution in the ISM 
with a large-scale density gradient. We do not involve into 
the analysis the small-scale emission structures which are 
seen on the SNR's maps, because these structures are re- 
sponsible for the local rises and do not essentially modify 
the global behaviour of the surface brightness. The plasma 
emission model used here is CIE. 



2. The 2-D Modelling of RCW86 

2.1. Method and models 

Evolution of the SNR in the medium with a large-scale 
density distribution is taken as the principal basis for the 
models presented in this paper. We accept here spherically 
symmetrical energy emission during a SN explosion. We 
assume that the nonspherical shape and the global regular 
surface brightness anisotropy of the SNR is caused only by 
the large-scale stucture of nonuniform medium. 

The new model parameters associated with orientation 
of the SNR and his projection on the plan of the sky ap- 
pear when 3D modelling executes. Projection effects com- 
plicate the analysis of the observations. In general case, it 
is impossible to reproduce unambiguily the real 3D shape 
for the knowing projection. In case of considered here the 
axis-symmetrical models, projection effects can easily be 
included into model. Simple connections between the pa- 
rameters of the model and the projection exists. The angle 
6 between the symmetry axis of the 2D SNR and the vi- 
sual cross-section is a new additional free paremeter of the 
model. 

At our previous paper [Paper 7) it have been described 
a new approximate analytical method which allows the 
modelling a point explosion in media with a large-scale 
density gradient. We apply this method here to the dis- 
cribing of an evolution of the SNR RCW86. The possi- 
bilities of the SN explosion in 185 A.D. and in the OB- 
assosiation will be considered separately . 

We take as the basic initial model's parameters the 
three observational results: the visual angular size of the 
SNR 8 « 40', the temperature of the plasma's hot com- 
ponent Thigh = 3.4 ± 0.2 keV, which corespond to cor- 
rected for the contribution of the galactic ridge emission 
spectrum from the whole SNR (Claas et al. 1989) and the 
surface brightness distribution profile al ong th e symmetry 
axis of the SNR's image (Pisarski et al. |l984[ ) (see Fig. |[ 
0). This profile is sensitive to large-scale density distribu- 
tion of ISM and allows to see the properties of the density 
gradient. 
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Table 1. Parameters of RCW86 derived from the X-ray observations. R is radius of the SNR, E^i is the energy of SN explosion 
in erg, is hydrogen number density of surrounding medium, M is swepted up mass, M© is the mass of Sun. 



Hydrodynamics model" 


PEM'' 


t, years 


R, pc 


d, kpc 


£51 




Wh, cm 


M, Mo 


Ref.'' 


Free expansion (SN II) 


CIE 


1800 


9 


1.4 


0.1 ^ 


0.2 


0.13-0.25 


7^ 14 


[1] 


Sedov (SN I) 


CIE 


1800 


6.5 


1.4 


0.1 -f 


0.2 


0.3 


5 


[1] 


Free expansion or Sedov 


NEI 


400 ^ 1900 


2.5 H- 9 


0.4^ 1.6 


0.1 ~ 


0.2 


0.04 H- 0.2 


0.1 ~ 21 


[2] 


Sedov for NE part of the SNR 


NEI 


1800 


6.7 


1.1 


0.15 




0.11"= 


4.9 


[3] 



^ All model assumed uniform ISM (density of ISM p° = const) 

^ Plasma emission model 
The hydrogen number density of the ISM for the SW part is es timated as {n^)sw 
Reference: [1] - Pisarski et al. (|l984[); [2] - Nugent et al. (|l984|); [3] - Claas et al. (|l989|). 



10(n^) 



These characteristics are supplemented additionally 
with the fourth one depending on the one from two basic 
assumptions: the first assumption is that RCW86 is a re- 
sult of the SN explosion in 185 A.D. (that gives the age of 
the SNR) and the second one is that RCW86 is a result 
of the SN explosion in the OB-assosiation (this yields the 
distance to the SNR). 

Th e temperature Tg at the shock front in the Sedov 
(|1959D model o f the S NR may be obtancd from Thigh: « 
Thigh/l-S (Itoh [l977 ). In Paper I it have been shown that 
the nonsperical SNR may be characterised by the some 
"average" characteristic temperature Teh oc {M)~^ , which 
will be likely connected with Thigh : 



Tch ~ Thigh/1.3 = (3.03 ± 0.2) • 10^ K. 



(1) 



We have also shown that during the adiabatic stage the 
value of Tch is close to Tg of the Sedov SNR with the same 
initial model parameters (e.g., the energy of explosion Eg, 
the hydrogen number density in the place of explosion 
71^(0)) with the maximal error about 20-;-30%. Therefore, 
we may write analoguously to the Sedov case that 



T 



ch 



T = 2.08 • 10 



11 



iH(0) 



2/5 



-6/5 

yr 



K, 



6.47 • lO-' 



E. 



51 \ /-r; 



(i?pc)"' K, 



(2) 



(3) 



where E^i — Eo ■ 10~^^, tyr is the age of the SNR in years, 
Rpc is the average radius of the nonspherical SNR in pc. 
Hereafter 7 = 5/3. 

For the calculation of the X-ray e missio n of the plasma 
in the CIE the Raymond & Smith ( 1977 ) data were ap- 
proximated. 



2.2. RCW86 as the result of the SN explosion in 185 A.D. 
2.2.1. Exponential medium 

Let us take as initial SNR characteristics the angular size 
6, the temperature Tch, the surface brightness profile and 
additionally the age of the SNR t = 1800 years. 



We will consider the evolution of the SNR in the 
nonuniform medium with the flat exponential density dis- 
tribution 



p°(r,^) =p°(O)exp(--^cos0) 



(4) 



where p°(0) is the initial density around the place of SN 
explosion, H is the scale of the density nonunifirmity, 9 is 
an angle between the maximal density decreasing direc- 
tion and the considered direction. Such exponential dis- 
tribution is the one of the most possible to approximate 
real continuous density distributions in interstellar clouds, 
gaseous galactic disc etc. 

When we fixed Tch and t we have got from @) an 
estimation on 



^5i/nH(0) = 1.5 ±0.25. 



(5) 



From this relation and appropriate explosion energy range 
i?5i — 0.1 1 the range for the initial number density 
n° (0) = 0.06 0.8 cm-3 follow. 

From the definition of a dimentionless time (Hnatyk & 
Petruk p^) T = i/im with = al"^ E'^'"^ pl'"^ {Q)rW^ 



(where aA is the self-similar constant, i?m is the distance 
scale) it is obtaned for 7 = 5/3 that 



t = rt,^ = 18.01 



^51 

«h(0) 



-1/2 



R 5/2,- vr 



(6) 



Rra,-pc is the distanse scale in pc. 

From this relation and (j^) a connection between the 
scale high H and the dimentionless time follows: 



H = i?,n =i (6.8 ± 0.2) • T 



'2/5 



pc. 



(7) 



The possible values of H lie in the range H = 150 pc 
for the Galactic's disk to few parsecs for a local nonunifor- 
mity of the ISM. The influence of the nonuniform medium 
on the dynamics of the SNR become more apprecieble with 
the age of a SNR. The surface brightness distribution of 
a SNR will be close to the Sedov one for small r. Since 
the observational distribution of the surface brightness in 
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Fig. 1. Surface brightness 5* distribution in the diapason 
e = 0.1 2 keV along of the median of the SNR in the ISM 
with the exponential density distribution for the four cases 
of the initial model parameters. It is suppose here for all cases 



that £;5i/n?[(0) = 1.5, K^^t 



2/5 



6.85, t = 1800 years and 



Teh = 3.0 -10 K. Solid lines show the cases with H = 70 pc, 
dashed - the cases with H — 10 pc. Energies of explosion; a, 
b: Esi = 1; c, d: E51 = 0.15. Dots represent the cases c, d 
muhiplied by the factor {{n'^{0))top/{n'^{0))tottomf = 10^ ^^ 



RCW86 is not similar to the Sedov one, the dimentionless 
time r have not to be small. 

The observed ratio ^ between the axes of the visual 
shape of RCW86 is in range £, = 1.0 1.2 (Pisarski et 
al. |1984 ). This corresponds to ranges of r = 0.01 ^ 11 
(Fig. 5 in Paper I). The scale hight H, as it follow from 
then must be anywhere between 2.6 and 40 pc. Such 
wide range of is a result of very slow decreasing of the 
visual shape's anizotropy with time (the same Fig.). Small 
scale hight (like H = 2.6 pc) it seems to be impossible be- 
cause the surface brightness contrast then must be about 
10^'^ times. Possible presence in the shape anizotropy the 
residual component related to anizotropy which the SNR 
might have on the free expansion stage (due to ashperical 
explosion) is the second reason which complicates using 
the small anizotropy of the shape for estimating the H 
and r. 

Physically the most correct limitation on H and t wc 
may get from the profiles of the surface brightness. In 
Fig. |l| are shown the profiles of the surface brightness for 
the four cases of the initial model peremeters. Relations 
(H) and (0) have place for each set of the parameters. We 
may see that for fixed H and r the shape of the profile 




10.0 15.0 
H, pc 



20.0 



Fig. 2. Ratios of the values of the two peaks in the surface 
brightness S distribution versus the high scale H for models 
of SNR in the flat exponential media (m. Here t — 1800 years 
and £51/71^(0) — 1.5. In the figure both contrast of surface 
brigh tness are shown: smoothed to 3' like observed (Pisarski 
et al. 1984) and without smoothing. 



is the same. Value of the blast energy E51 (and as conse- 
quence of (^), the number density n^{0)) give the ampli- 
tude of the profile which may differ with the factor about 
few ten times. It is follow from Fig. |^ that the values H 
and T affect on the contrast of the surface brightness. 

For finding H, it is possible to accept the contrast be- 
tween the brightness in the different points of the SNR's 
image. For example, it may be used the relation between 
the two "near- front" maximums in the distribution. The 
most appropriate value of high scale for this choice, as 
it may be seen from Fig. ^, is = 11.0 pc (and respec- 
tively T — 0.3). Then calculations show (Fig. ^) that the 
most appropriate value of the blast energy is E51 = 0.22 
(7i;^(0) 0.15 cm-3). 

For normalizing the surface brightness S in photon en- 
ergy range e = 0.1-^2.0 keV, the effective values of photon 
energy e = 0.1 keV and absorption cross section a as cross 
section at photon energy e — 0.316 keV from Morrison & 



McCammon (1983) have been used. 

Obtaned SNR's luminocity i^'^^^ = 3.5 • lO^^ erg/s 
is consistent with esteblished by Pisarski et al. (1984) 
L^-2-4^2.1034<,erg/s. 

In Fig. H it is shown the maps of the surface brightness 
distribution for this model in comparison with the distri- 
bution of the spectral index a = —d\nF^/d\ne, where 
is the continuum flux at photon energy e. 

2.2.2. Power law medium 

It is useful to note, that concreete appearance of surround- 
ing density distribution is not strongly essential then size 




Fig. 3. The same as in Fig. ^ for tlie ratios between tlie max- 
imum in the surface brightness distribution Smax,i and the 
value of the brightness in the visual geometric center of the 
SNR Sc. 



of SNR does not exceed a few scale hights. Namely, similar 
parameters of model may be obtained under assumption 
that the SNR evolve in other type of media, if contrast of 
density along the surface of the SNR will be similar. Let 
us consider, for example, a medium with the spherically- 
symmetrical power law density distribution created by 
steller wind 



(8) 



then the SN explosion position is displaced on the distance 
To from the center of symmetry f — 0, therefore, the den- 
sity distribution as a function of the distance r from the 
explosion point r = is: 



p°(r,0)=p°(O)| 



2rro cos I 



(9) 



We take uj = —2 and calculate a number of models in 
order to get an appropriate using the surface brightness 
contrast between the two maximums in the brightness dis- 
tribution. 

Observed surface brightness contrast is obtaned in this 
model with = 21.8 pc. Characteristic scale Hch for 
medium @ with such Tq may be estimated. It equals 
Hch = 11 pc. Other parameters of the model are really 
close to the same parameters of the SNR in exponential 
medium (Q) and shown in Tables || and Fig. ^. 

2.2.3. Exponential plus uniform medium 

RCW86 may evolve in the ISM with the contact between 
the uniform medium p° = const and the region with the 



Fig. 4. Observed (line 1) distribution of the surface bright- 
ness S along of the symmetry axis of RCW86 (Pisarski et al. 



1984) and the distributions in the models of SNRs with the 
age t = 1800 years: line 2 - the exponential medium (^) with 
H = 11 pc and Esi = 0.22, n^{0) = 0.15 cm"^; line 3 - the 
power law medium M with To = 21.8 pc and the same £51, 



71^(0); line 4 - the exponential medium (^) with iJ" = 5 pc and 
£51 = 0.17, n^{0) — 0.11 cm~'^; line 5 - the uniform medium, 
£51 — 0.17, = 0.11 cm^'^. All distributions are smoothed 
to 3'. 



higher density (e.g., interstellar cloud) where density is 
distributed according to the exponential law (|^). In this 
case one part of the SNR (the NE part of RCW86) will 
evolve in the uniform medium whereas the evolution of the 
another part (the SW part of RCW86) will be determined 
by the variation of the density in the dense region. In 
this composite model the surface brightness in the central 
part of the SNR will be close to the brightness which the 
Sedov model gives, but the maximal brightness will be 
determined by the concrete distribution of the density in 
the nonuniform region. 

When we take as initial parameter the ratio of the 
maximum of the brightness distribution to the value of 
the brightness in the visual geometrical centre of the SNR, 
we obtaine another model parameters for the case of the 
exponential law density distribution, namely, H = 5.0 pc 
(Fig. |) and E51 = 0.17 (Fig. §. We may calculate now 
the profile of the surface brightness for the Sedov model 
wich discribes the NE part of RCW86 in the assumption 
of the composite medium (Fig. ^). The initial parame- 
ters for it are the same as for the SW part, E^i = 0.17, 
n^(0) = 0.11 cm^'^. We may see that the composite 
medium model (exponential plus constant density dis- 
tributions) more better describes the observational sur- 
face brightness distribution than the exponential medium 
alone. 
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Fig. 5. a The surface brightness S (in erg s^^ cm^^ st^^) dis- 
tribution in the photon energy range e = 0.1 — 2 keV. b The 
surface distribution of the spectral index q at e = 5 keV for the 
SNR model in the exponential medium (^) with _ff = 11.0 pc. 
Parameters of the model are t — 1800 years, E^i = 0.22, 
n^(0) = 0.15 cm"^ and 6 = 0°. 



It is nesessary to note, when RCW86 is the SNR cre- 
ated by the explosion in 185 A.D., we see it near the max- 
imum disclosing projection {6 « 0°). Really, another pro- 
jections hide the real contrasts (e.g., of axes size, surface 
brightness; see Paper I). If we see the contrast of surface 
brightness decreased already, then the real one must be 
greater and therefore H must be smaller. 

Obtaned parameters of RCW86 are summarised in Ta- 
ble H (the first three columns). 

2.3. RCW86 as the result of the SN explosion in the OB- 
association 

Let us suppose now that RCW86 is a result of the SN 
explosion in the OB-assosiation which distance from us is 



Table 2. Parameters of RCW86 obtaned from the different 
models of the SNR. Plasma emission model is CIE. Lx is in 
the photon energy range e = 0.1 -^2.0 keV. R{0) and ii(7r) are 
the radii of shock in directions ^ = and 9 = n respectively, 
D is the shock velocity. 



Parameters 




Initial density distribution 




of model 




PL" 


EU'' 






Observed 












O, arcmin 


40 


40 


40 


40 


40 


Teh, 10' K 


3.0 


3.0 


3.0 


3.0 


3.0 


Additionally supposed 










t, years 


1800 


1800 


1800 


— 




d, kpc 


— 


— 


— 


2.8 


2.8 


H, pc 


11 


~ IL* 


5.0 


26 


22 


5, degree 














30 


Obtained 












t, years 


1800 


1800 


1800 


4300 


4300 


So, 10^1 erg 


0.22 


0.22 


0.14 


2.0 


2.0 


n?[(0), cm"^ 


0.15 


0.15 


0.10 


0.1 


0.1 


ii(0), pc 


7.50 


7.42 


6.84 


17.96 


18.25 


R{-k), pc 


6.27 


6.20 


5.78 


14.96 


14.78 


D(0), km/s 


1800 


1740 


1490 


1810 


1870 


D{n), km/s 


1260 


1220 


1070 


1250 


1230 


r(o), 10' K 


4.45 


4.18 


3.03 


4.50 


4.82 


T(7r), 10' K 


2.17 


2.03 


1.58 


2.16 


2.07 


n'^{R,Q), cm"^ 


0.08 


0.08 


0.10 


0.18 


0.19 


nH(-R,7r), cm"^ 


0.27 


0.29 


0.31 


0.05 


0.04 


d, kpc 


1.18 


1.17 


1.08 


2.83 


2.83 


M, M0 


6.9 


6.8 




62 


62 


Zp(Lx, erg/s) 


34.54 


34.53 




35.33 


35.35 


a(5 keV) 


1.90 


1.88 




1.91 


1.92 



^ E is the flat exponential medium (|^ 
^ PL is the power law medium (^ with uo = —2 and 
ro = 21.8 pc 

EU is the exponential medium (El) plus uniform 



Characteristic scale hight _ffch 



estimated (e.g., Westerlund 1969). In this section the next 
question is the main: what scale hight has the medium 
around the RCW86 if the SNR's progenitor have exploded 
in this OB-assosiation? Therefore the nonuniform media 
with exponential density distribution (^ is used. 

We take as the initial SNR characteristics the obser- 
vational angular size 6, the temperature Teh) the surface 
brightness profile and, additionally, the distance to the 
remnant d w 2.8 ± 0.4 kpc (Rosado et al. |l99§). 

It is possible now to estimate the average SNR's radius 
from e and d: Rp^ = 0.5dkpce73.438 w 16.3 ± 2.3 pc. 
Therefore we have got from (|^) an estimation on 



E5i/n'^{0) ~ 23 ± 10. 



(10) 



This means that the explosion energy must be high and 
the initial number density have to be low. 

The age of the SNR as it follows from (||) is any- 
where from t = 3500 years for Teh — 2.8 • lO' K and 




Fig. 6. Ratios of tiie values of tiie two peaks in tlie X-ray 
surface brigiitness S distribution versus H for the SNRs in the 
exponential media when t = 4300 years, ^51/72^(0) = 20.3 
and 5 = 0°. Dashed line is the same relation for the same SNR 
models, except S = 30°. All ratios are smoothed to 3'. 



Fig. 7. Distribution of the surface brightness S along the sym- 
metry axis of the SNR RCW86: line 1 - observational (Pisarski 
et al. 1984); line 2 - model in exponential medium with 
H = 25.7 pc and t = 4300 years, E^i = 2.0, n'^{0) = 0.1 cm"^, 
S — 0°; line 3 - the same model, except H — 22.3 pc and 
(5 = 30°. All distributions are smoothed to 3'. 



E5i/n^{0) = 13, up to t = 5350 years in case Teh = 
3.2 ■ 10'' K and £^5i/ng,(0 ) = 3 2. We may see, with agree- 
ment with Rosado et al. ( 1996|) , if RCW86 is the remnant 
of the SN exploded in OB-assosiation it can not be borned 
in 185 A.D. 

For the distance d = 2.8 kpc and Teh = 3.0 ■ 10^ K 
the age t — 4300 years and ratio E5i/n^{Q) ~ 20.3 have 
been obtained. We calculate the evolution of RCW86 in 
the exponential medium (Q) with these initial parameters. 

It may be seen from Fig. ^ that the scale hight H — 
25.7 pc is the most appropriate for the case of the maximal 
disclosure of the remnant {6 — 0°). 

In order to estimate an explosion energy Eo and a num- 
ber density n^{0) it was calculated a number of the models 
with such H and different Eq. It may be seen from Fig. [t] 
that the surface brightness distribution along the symme- 
try axis for £'51 = 2.0, n^{0) =0.1 cm~^ is near of the 
observed one. The whole maps of the surface brightness S 
and the spectral index a distribution are shown in Fig. |[ 

Fig. ^ demonstrates also the surface brightness distri- 
bution along the symmetry axis for the case of the dif- 
ferent projection of the SNR, when S — 30°. We see that 
distribution, like observational, may be obtaned also un- 
der different projection of the 2D SNR on the plan of the 
sky with the smaller H = 22.3 pc. Now we cannot separate 
one case of projection from another but we see also that 
the model with S — 30° does not essentially change the 
SNR's parameters which are summarised for these differ- 
ent models of RCW86 in the last two columns of Table |[ 

Modelling reveal also that for model with 6 = 30° both 
the surface brightness distribution and the surface distri- 



bution of the spectral index are also very close to shown 
in Fig. H case of (5 = 0°. 

2.4- Contribution of the ejecta 



Nugent et al. (1984) have made the conclusion that the 
observed data on RCW86 can be fitted without taking into 
consideration the emission from the ejecta. These authors 
have also suggested that the value of f3 — M^-J E^i (Mej 
is the mass of ejecta in Mq) lies much probably in the 
range 1 to 10. In our models of the SNR the blast energy 
are £'51 = 0.22 in the case of SN explosion in A.D. 185 
and £51 = 2.0 in the case of SN explosion at a distance 
d = 2.8 kpc. Terefore, the mass of ejecta may be estimated 
as Mcj = 0.22 2.2M0 or M^j = 2 ^ 2OM0. Swepted up 
masses are respectively 7 and 62 M© and exceeds A/cj in 3 
to 31 times. We expect that under our assumption of CIE 
conditions and obtained ~ 0.1 cm^'^, such swepted 
up masses are enough for RCW86 to be in Sedov phase 
of his evolution. So, we may assume that in both cases 
the ejecta do not considerably modify the emission from 
RCW86. However, if to take into account the NEI effects 
and consider the value of the (3 up to the possible upper 
limit which have been put as [3 ^ AQ, the contribution of 
the ejecta into the emission might be essential (Nugent et 
al. |1984D . On the other hand, if RCW86 would evolve in 
10^), it would entry in Scdov stage 



very dense ISM (n^ ^ 
after swepting up M 
|T996| ). 
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It is shown that observed surface brightness distribu- 
tion of RCW86 may be obtaned in the models with the 
two different initial assumptions: one about the Supernova 
explosion in 185 A.D. and another about the explosion in 
the OB-asossiation. Data we posses do not allow us to 
diside which of those models is true. The parameters ob- 
taned for RCW86, basing on these two assumption, are 
summarased in Table ^. The models give the observational 
contrasts of the surface brightness when to consider the 
ISM with the scales of nonuniformity 11 pc (if the age of 
RCW86 is 1800 years) or 20 - 25 pc (when the SNR is dis- 
tant from us on 2.8 kpc). The preshock density contrasts 
between the southwestern and the northeastern parts of 
RCW86 in range 3.5 — 4.5 are obtained for both considered 
assumptions. 
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Fig. 8. a The surface brightness S (in erg s~^ cm~^ st~^) dis- 
tribution in the range 0.1 — 2 keV and b the surface distribu- 
tion of the spectral index a at e = 5 keV for the SNR model in 
the nonuniform exponential medium (^ with H = 25.7 pc. 
Parameters of the model are t = 4300 years, £51 = 2.0, 
n?[(0) = 0.1 cm-3 and S = 0°. 



3. Conclusions 

Observations reveal complicated RCW86 morphology. 
Close to spherical shape of the SNR coexists with the sur- 
face brightness distribution which is far from the Sedov 
one. When we take into account the nonuniform ISM we 
may explain such morphology. 

We considered here the 2-D models of RCW86 in a few 
types of nonuniform media. It is shown that the features 
of the surface brightness distribution allow us to restore 
the SNR characteristics. As it is described, the contrasts 
in the surface brightness depend on the gradient of the 
density of the ISM and the age of the SNR, thereas the 
amplitude of the surface brightness distribution depends 
on the energy of the Supernova explosion and the initial 
density around the place of explosion. 
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